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Balance functions have been measured in terms of relative pseudorapidity (A rf) for charged particle 
pairs at the Relativistic Heavy-Ion Collider (RHIC) from Au+Au collisions at ^/snn = 7.7 GeV to 
200 GeV using the STAR detector. These results are compared with balance functions measured 
at the Large Hadron Collider (LHC) from Pb+Pb collisions at ^/snn = 2.76 TeV by the ALICE 
Collaboration. The width of the balance function decreases as the collisions become more central 
and as the beam energy is increased. In contrast, the widths of the balance functions calculated 
using shuffled events show little dependence on centrality or beam energy and are larger than the 
observed widths. Balance function widths calculated using events generated by UrQMD are wider 
than the measured widths in central collisions and show little centrality dependence. The measured 
widths of the balance functions in central collisions are consistent with the delayed hadronization of 
a deconfined quark gluon plasma (QGP). The narrowing of the balance function in central collisions 
at ^/snn = 7.7 GeV implies that a QGP is still being created at this relatively low energy. 


PACS numbers: 25.75.Gz 
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Event-by-event charge correlations and fluctuations 
can be used as a tool to study the dynamics of hadroniza- 
tion in relativistic heavy-ion collisions DJ-l33j. One such 
observable, the balance function |27 h 30I ]. is sensitive to 
the correlation of balancing charges. The basic idea of 
the balance function is that charge is created in balanc¬ 
ing pairs that originate from the same point in space and 
time. By means of a like-sign subtraction, the balance 
function can yield the distribution of relative momen¬ 
tum between the balancing charges. Balance functions 
are sensitive to the mechanisms of charge formation and 
the subsequent relative diffusion of the balancing charges 
f+zt and are also affected by the freeze-out temperature 
and radial flow [28[. Model calculations show that collec¬ 
tive flow is not sufficient to explain the balance-function 
widths measured in central Au+Au collisions at ^/snn 
= 200 GeV [H, (33, (HJ. Balance functions for central 
collisions have been shown to be consistent with blast- 
wave models where the balancing charges are required to 
come from regions with similar collective flow [30J- The 
inferred high degree of correlation in coordinate space_has 
been postulated as a signal for delayed hadronization [27j ■ 
In central collisions, a deconfined system of quarks and 


gluons is created, which cools and expands [35]. Most 


of the observed balancing charges are then created when 
the deconfined system hadronizes, which limits the time 
available for the balancing charges to diffuse away from 
one another. This leads to tighter correlations in coor¬ 
dinate space of balancing charges, and due to collective 
motion, results in tighter correlations in relative momen¬ 
tum and relative rapidity. Alternatively, if the charges 
are created early (on the order of 1 fm/c), the balancing 
charges are less correlated in the final state because the 
balancing charges have more time to move apart from 
one another. Thus, a narrow balance function in terms 
of relative pseudorapidity or relative rapidity in central 
collisions compared with peripheral collisions implies de¬ 
layed hadronization. 

The balance function is a conditional distribution f2?| , 
which can be written as: 

| N-+(Ar)) - AT—(At?) \ ^ 


The balance function in terms of Ar) (B (A?;)) represents 
the probability of seeing a particle which has a relative 
pseudorapidity A?y with respect to its opposite sign part¬ 
ner, given the condition that its opposite sign partner 
has already been seen inside the detector. Specifically, 
Ak|_(Ary) is the histogram of Ary for all the positive¬ 

negative charged particle pairs in an event. It is then 
accumulated over all events. Similarly, the histograms 
for V ++ , N _, and N-+ are calculated. For the de¬ 

nominators, N .is the number of positive (negative) 
particles integrated over all events. 

The system size and centrality dependence of the bal¬ 
ance function for all charged particles has been studied 
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2010 

32 
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28 
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2011 
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2010 

7.7 

7.7 

2010 

2.2 


TABLE I. Summary of the data used in this analysis. 


by the NA49 Collaboration at ^/snn = 17-3 GeV for 
p+p, C+C, Si+Si, and Pb+Pb collisions [36j- The bal¬ 
ance function for all charged particles narrows in cen¬ 
tral Pb+Pb collisions at 17.3 GeV and the widths of 
the balance functions for p+p, C+C, Si+Si, and Pb+Pb 
collisions scale with the number of participating nucle¬ 
ons. The NA49 Collaboration has also published results 


37[ for the rapidity dependence and beam energy depen¬ 


dence of the balance function for all charged particles for 
Pb+Pb collisions from ^/snn = 6.3 GeV to 17.3 GeV. 
The balance function was observed to narrow in central 
collisions for midrapidity, but did not narrow at forward 
rapidity. The authors of Ref. ]37f showed that the nar¬ 
rowing of the balance function in terms of Ary in central 
collisions was explained with the AMPT (a multiphase 
transport) model ]38) incorporating delayed hadroniza¬ 
tion, while models such as HIJING (heavy-ion jet inter¬ 
action generator, version 1.38, default parameters) Q 
and UrQMD (ultra relativistic quantum molecular dy¬ 
namics, version 3.3, with default parameters) [40] failed 
to reproduce the observed narrowing. 

The STAR Collaboration has presented a study of the 
longitudinal scaling of the balance function in Au+Au 
collisions at ^/snn = 200 GeV |4lj |. STAR has pub¬ 
lished results for balance functions from p+p, d+Au, 
and Au+Au collisions at ^/snn = 130 and 200 GeV in 
terms of Ary, relative rapidity (Ary), relative azimuthal 
angle (Ac/>), and invariant relative momentum (q; nv ) for 
all charged particles, for charged pions, and for charged 
kaons [42[ S3 ■ The balance functions for all charged par¬ 
ticles and for charged pions narrow as the events become 
more central while balance functions calculated using 
HIJING and UrQMD showed no centrality dependence. 
The ALICE Collaboration has recently published mea¬ 
surements [HJ from Pb+Pb collisions at ^/snn = 2.76 
TeV that also show that the balance functions in terms 
of Ary narrow in central collisions. 

In this paper, we report measurements of balance func¬ 
tions for all charged particles with 0.2 < px < 2.0 GeV/c 
in terms of relative pseudorapidity (Ary) from Au+Au 
collisions at ^/snn = 7.7, 11.5, 19.6, 27, 39, 62.4, and 
200 GeV. We observe that the balance functions narrow 
in central collisions and narrow as the beam energy is 
increased. We compare with the results from ALICE 
for Pb+Pb collisions at ^/snn = 2.76 TeV by restricting 
STAR’S acceptance to Ary < 1.6 to match the acceptance 
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FIG. 1. (Color online) The balance function in terms of A?y 
for all charged particles with 0.2 < pt < 2.0 GeV/c from cen¬ 
tral Au+Au collisions (0-5%) for ^/snn from 7.7 to 200 GeV. 
The data are the measured balance functions corrected by 
subtracting balance functions calculated using mixed events. 
Also shown are balance functions calculated using shuffled 
events. 


of ALICE, correcting for the acceptance of STAR in rj , 
and calculating the width of the balance function over the 
range 0.1 < Ary < 1.6, where the lower limit of 0.1 is cho¬ 
sen to suppress effects from interpair correlations [e.g., 
Hanbury-Brown and Twiss (HBT) and final-state inter¬ 
actions]. We observe that the balance function in terms 
of Ary narrows as the beam energy is raised from ^/snn = 
7.7 GeV to 2.76 TeV. When the observed balance func¬ 
tion widths are scaled by the width observed in the most 
peripheral bin, the relative widths still decrease as the 
events become more central and as the beam energy is 
increased. These results contrast with those presented 
in Ref. (44[, where the scaled balance function widths 
are shown to be nearly the same at RHIC and LHC en¬ 
ergies. The present observations are consistent with the 
concept of delayed hadronization of a deconfined quark 
gluon plasma (QGP) with the deconfined system having 
a longer lifetime at the highest energy. The narrowing of 
the balance function in central collisions at ^/snn = 7.7 
GeV implies that a QGP might still be created at this 
relatively low energy. 

The data were taken with the STAR detector [45]] dur¬ 
ing the years 2010 and 2011. TableUshows a summary of 
the data sets used in this analysis. Au+Au collisions were 
studied at seven beam energies ranging from 7.7 GeV to 
200 GeV. The centrality of each collision was determined 
according to the measured charged hadron multiplicity 
within the pseudorapidity range |ry| < 0.5. Nine central¬ 
ity bins were used: 0-5% (most central), 5-10%, 10-20%, 
20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80% 
(most peripheral). At each of the seven beam energies, 
the average number of participating nucleons, V part , is 
calculated for each of the nine centrality bins using a 
Glauber model. To ensure a more uniform detector ac¬ 
ceptance, events were accepted only when the position 
of the reconstructed primary vertex was within 30 cm of 
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FIG. 2. (Color online) Energy dependence of the balance 
function widths compared with the widths of the balance 
functions calculated using shuffled events. Also shown are 
the balance function widths calculated using UrQMD. The 
dashed line represents the width of the balance function cal¬ 
culated using shuffled events for a constant dN/di 7 distribu¬ 
tion. Error bars represent the statistical error and the shaded 
bands represent the systematic error. 
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FIG. 3. (Color online) Balance function widths for the most 
central events (0-5%) compared with balance function widths 
calculated using shuffled events. Also shown are balance func¬ 
tion widths calculated using UrQMD and shuffled UrQMD 
events. The dashed line represents the width of the balance 
function calculated using shuffled events for a constant dN/di 7 
distribution. 
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the center of STAR (|Vertex| < 30 cm). In addition, the 
radial position of the primary vertex was required to be 
less than 2 cm from the center of the beam line to avoid 
beam pipe events. All events were required to have at 
least one matched track with the STAR Time-of-Flight 
(TOF) system [|(| to suppress pile-up events. 

All tracks in the Time Projection Chamber (TPC) 
were required to have more than 15 measured space 
points along the trajectory. The ratio of the number 
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FIG. 4. (Color online) Acceptance-corrected balance function 
widths for Au+Au measured over the range 0.1 < A?; < 1.6 
compared with similar results from Pb+Pb collisions from 
ALICE [+|. Only statistical errors are shown. Lines repre¬ 
sent fits of the form a + b(N pa , r t) 0 ' 01 . 


of reconstructed space points to possible space points 
along the track was required to be greater than 0.52 to 
avoid track splitting. Tracks in the TPC were character¬ 
ized by the distance of closest approach (DCA), which is 
the smallest distance between the projection of the track 
and the measured event vertex. To suppress decay ef¬ 
fects and background, all tracks were required to have 
a DCA less than 3 cm. A transverse momentum cut 
of 0.2 < pt < 2.0 GeV/c and a pseudorapidity cut of 
1771 < 1.0 were applied. 

In addition to real data, mixed events and shuffled 
events were also used in this analysis. Mixed events are 
created by grouping the events according to bins in cen¬ 
trality and bins in the position of the reconstructed ver¬ 
tex of the event along the beam direction. Ten centrality 
bins and five bins in z ver t ex were used. A set of mixed 
events is created by taking one track chosen at random 
from an event, which is selected according to the bin in 
centrality and the bin in event vertex position. A mixed 
event includes no more than one track from any observed 
event. This mixed-event data set has the same number 
of events with the same multiplicity distribution as the 
original data set but all correlations are removed. The 
mixed-event subtraction was important, especially at low 
energies, to account for the effects caused by unbalanced 
positive charges in each event. 

Shuffled events are produced by randomly shuffling the 
charges of the particles in each event, which removes the 
charge correlations while retaining global charge conser¬ 
vation. Because shuffling uniformly distributes a parti¬ 
cle’s balancing partner across the measured phase space, 
balance functions calculated using shuffled events can be 
used to gauge the widest balance functions that one can 
measure within the experimental acceptance of STAR. 

Fig. |T| shows the balance functions in terms of Ap for 


all charged particles for Au+Au collisions at ^/snn = 7.7, 
11.5, 19.6, 27, 39, 62.4, and 200 GeV for the most central 
events (0-5%) along with balance functions calculated us¬ 
ing shuffled events and balance functions calculated using 
mixed events. The data shown in the figure are the mea¬ 
sured balance functions corrected by subtracting balance 
functions calculated using mixed events. These data have 
not been corrected for efficiency or acceptance. The con¬ 
clusions of this paper involve the width of the balance 
function in which the efficiency cancels out. The model 
calculations shown in this paper use the STAR accep¬ 
tance. When comparisons are made withjthe width of 
the balance functions reported by ALICE [44| , the STAR 
data are corrected for acceptance. 

At the lower energies, the balance functions calculated 
using mixed events exhibit an oscillatory distribution 
that is 0 at A 77 = 0, has a positive value at Ap = 0.5, is 
0 at Ap = 1, has a negative value at A 77 = 1.5 and is 0 
again at Ap = 2. This oscillatory behavior lessens as the 
events become more peripheral and as the beam energy 
is increased. This effect is due to unbalanced positive 
charge that is not subtracted by the same sign subtrac¬ 
tion inherent in the balance function. The additional 
positive charges are dominantly protons and have a dif¬ 
ferent dN/dp distribution than the negative charges that 
are dominantly pions. The dN/dp distributions for the 
difference between the positive and negative charges have 
minima at ?y = -1, p = 0, and p = 1. Thus, when the 
balance function in terms of A p is calculated for mixed 
events at the lower energies and in more central colli¬ 
sions, the oscillatory distribution is obtained. At ^/snn 
= 200 GeV, the balance functions calculated using mixed 
events are zero for all centralities, which indicates that 
the amount of unbalanced positive charge is small. As 
the beam energy is decreased, the unbalanced positive 
charge increases and the balance functions calculated us¬ 
ing mixed events become significant. 

The corrected balance functions are narrower than the 
balance functions calculated using shuffled events and 
the balance functions narrow as the events become more 
central (see below). Also visible are the effects of in¬ 
terpair correlations [HBT and final-state interactions] 
that model calculations have shown to be significant for 
A ?7 1$ 0.1 j29j. Specifically, B (Ap) for Ap < 0.1 is not- 
icably higher than the trend of the remaining points at 
7.7 GeV while B (Ap) for (A 77 ) < 0.1 is lower than the 
trend at 200 GeV. The width of the balance function is 
characterized in terms of a weighted average: 


(Ap) 


Supper 

£ B(Api)Api 

^ = *lower 

^upper 

E B(Api) 

X —Slower 


( 2 ) 


Here i is the bin number and B(Api) is the value of 
the balance function for the relative pseudorapidity bin 
Api. The weighted average is calculated over a range 
in Ap chosen to minimize contributions from HBT and 
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Coulomb effects (Ary >0.1) and maximize the acceptance 
of STAR (Ary < 2.0). 

Fig. [2] shows the balance function widths for Au+Au 
collisions from ^/snn = 7.7 GeV to 200 GeV for nine cen¬ 
trality bins. The widths are calculated for 0.1 < Ary < 
2.0 to remove the distortions caused by interpair corre¬ 
lations for Ary < 0.1 [29|. The widths of the balance 
functions calculated using shuffled events are larger than 
the widths of the balance functions calculated using data. 
The widths of the balance functions using shuffled events 
shown in Fig. [2] are close to the value 0.733, which one 
would expect for shuffled events from a flat dN/dr/ dis¬ 
tribution over the range — 1 < ry < 1. The data show a 
smoothly decreasing width with increasing beam energy 
and as the collisions become more central. Fig. [5] also 
shows the widths of balance functions calculated using 
UrQMD. The UrQMD calculations are analyzed in the 
same way as the data with the balance functions calcu¬ 
lated from mixed UrQMD events being subtracted from 
the balance functions calculated using UrQMD. For beam 
energies below 20 GeV, the balance function widths from 
UrQMD increase as events become more central whereas 
the measured widths decrease. Above 20 GeV, the bal¬ 
ance function widths from UrQMD show little centrality 
dependence. In peripheral collisions, the balance func¬ 
tions widths from UrQMD approach the value of the 
measured balance function widths. The UrQMD model 
is a hadronic model that does not have a deconfined 
phase and has little flow. The early hadronization time of 
the particles calculated using UrQMD combined with the 
strong interaction between final state particles causes the 
larger balance function widths in central collisions while 
the balance function widths calculated using UrQMD are 
close to the measured balance function widths in periph¬ 
eral collisions. 

One source of systematic errors was estimated by 
studying the difference between the 200 GeV results 
from three different runs (in 2007, 2010, and 2011) that 
used different tracking software and incorporated differ¬ 
ent hardware configurations in STAR. A second source 
of systematic errors was estimated by varying the DCA 
used to select tracks. A third source of systematic er¬ 
ror was estimated by varying the range of the Vertex of 
events accepted in STAR. The systematic errors in the 
extracted widths are shown as a shaded band in Fig. [2] 
Note that the systematic error in the width for the most 
central bin at all energies was of the same order or less 
than the statistical errors. 

Fig.[3]shows the width of the balance function in terms 
of Ary for central collisions (0-5%) as a function of beam 
energy. The measured balance function widths decrease 
smoothly with increasing beam energy. Also shown are 
the widths of the balance function calculated using events 
generated with the UrQMD model. Although the energy 
trends for the width of the balance function in UrQMD 
and data appear similar, the data are much narrower, 
and as shown in Fig. 2, UrQMD predicts the wrong cen¬ 
trality dependence. The widths of the balance function 



FIG. 5. (Color online) Acceptance-corrected balance function 
widths for Au+Au measured over the range 0.1 < Ary < 1.6 
normalized to the most peripheral centrality bin compared 
with similar results from Pb+Pb collisions from ALICE [44| . 
Only statistical errors are shown. Lines represent fits of the 
form a + b(N pa , lt ) 0 01 . 


calculated from shuffled events from both the data and 
UrQMD are much larger than the widths calculated us¬ 
ing the data. The decrease of the shuffled widths at the 
lower beam energies reflects the fact that the dN/drj dis¬ 
tributions at the lower beam energies are not completely 
flat. The fact that the measured balance function widths 
decrease smoothly with increasing beam energy and are 
much smaller than the widths predicted by UrQMD is 
consistent with the idea of delayed hadronization. 

To compare with the balance functions measured by 
ALICE, we correct our measured balance functions for 
the acceptance of STAR in ry using the expression [28j : 

7 ?Ar; max =2 (Ary) = Boo (Ary) ^1 , (3) 


where B oo (Ary) is the STAR balance function corrected 
for acceptance in ry assuming that STAR’s acceptance 
is constant for —1 < ry < 1, and i?A? 7 ma . x =2 is the 
measured STAR balance function that is not corrected 
for acceptance in ry. For the comparison with the re¬ 
sults from ALICE |3|, we calculate the widths of the 
acceptance-corrected balance functions over the range 
0.1 < Ary < 1.6 to suppress effects from interpair cor¬ 
relations and to match the acceptance of ALICE in ry. 
Fig. @] shows these widths as a function of centrality and 
beam energy for Au+Au collisions. In the same figure 
we show the width of the balance function from Pb+Pb 
collisions at ,/snn = 2.76 GeV calculated for three pub¬ 
lished centralities [H} over the same range in Ary. The 
balance functions narrow as the beam energy is raised 
from yj snn =7.7 GeV up to 2.76 TeV and the balance 
functions narrow as the collisions become more central. 
These observations are consistent with the concept of de¬ 
layed hadronization. 

The authors of Ref. 44] assert that the relative de- 
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crease of A?; with centrality does not change appreciably 
with beam energy. To address this point, we calculate 
the ratio of the width of the balance function at each cen¬ 
trality to the width of the balance function in the most 
peripheral bin at each beam energy, (A 77 } /(A? 7 } peripheral . 
Because the peripheral bin at the lower energies has low 
statistics, we first fit the measured widths at each beam 
energy with a function of the form a + b(N paTt ) 001 and 
then take the ratio of the measured widths to the width 
of the fitted distribution at the most peripheral centrality 
bin. These results are shown in Fig. [5j 

The relative decrease of the balance function width is 
much larger at ^/snn = 2.76 TeV. The relative decrease 
then gets smaller as the beam energy is lowered. Thus, we 
observe that the relative decrease of the balance function 
width clearly changes with beam energy. The difference 
between the present analysis and the one presented by the 
authors of Ref. [44[ is that we calculate the widths over 
the range 0.1 < A 77 < 1.6 for both experiments, which 
minimizes the contributions from interpair correlations 
and maximizes the measured acceptance. In contrast, 
the authors of Ref. [Hj calculated the widths over the 
range 0.0 < A 77 < 1.6 for the ALICE balance functions 
and 0.0 < Arj < 2.0 for the STAR balance functions. 
We do not compare with the results from NA49 [13] here 
because the acceptance of NA49 in 77 is relatively small. 

Model calculations [H, HH show that a part of the 
narrowing of the balance function in central collisions 
is due to radial flow. Thus, the fact that the balance 
functions in central Pb+Pb collisions at 2.76 TeV are 
narrower than those in central Au+Au collisions at 200 


GeV may be due to an increase in radial flow. One would 
expect that the balance function in terms of A 77 would 
be narrower for a longer-lived deconfined QGP, which 
implies that these results are consistent with the concept 
of delayed hadronization. 

In conclusion, we observe that the balance function in 
terms of A?/ is narrow in central collisions of Au+Au. At 
higher beam energies, the balance function in terms of 
A77 in central collisions of Pb+Pb is even narrower. This 
observed narrowing is consistent with the concept of the 
delayed hadronization of a deconfined QGP produced in 
these collisions. We observe that the balance functions 
in Au+Au events at ^/snn = 7.7 GeV still narrow as 
the collisions become more central, which suggests that a 
deconfined QGP might still be produced at this relatively 
low beam energy. 
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